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ABSTRACT 

Resultants are important special functions used in description of non-linear phenomena. Resultant Rri,....r„ defines a 
condition of solvability for a system of n homogeneous polynomials of degrees ri, . . . ,rn in n variables, just in the same way 
as determinant does for a system of linear equations. Unfortunately, there is a lack of convenient formulas for resultants when 
the number of variables is large. In this paper we use Cauchy contour integrals to obtain a polynomial formula for resultants, 
which is expected to be useful in applications. 



It is well-known, that a system of n linear equations in n variables has a non-zero solution if and only 
r^ , if the corresponding determinant vanishes. Solutions to non-homogeneous linear systems and Gaussian 

integrals (associated with linear equations of motion) are then expressed through determinants. Similarly 
in non-linear case a system of homogeneous polynomials /i(a?), . . . , fn{x) of degrees ri, . . . , r^ in n variables 
i-G ' X = {xi, . . . , Xn) has a non-zero common root, if and only if coefficients of the system satisfy a constraint 



-Rri,...,r„{/l, ■ • ■ ,/n} — 



where expression R is called a resultant of this system (the subscripts ri , . . . , r„ will be often omitted in 

what follows). For linear systems resultant is just a determinant, the well-known and well-understood object 

^ \ of linear algebra. Unfortunately, the knowledge and understanding of general resultants is not equally deep, 

yet. Their study is a relatively new subject, which is natural to call non-linear algebra [HIl]. 
\j^ \ In analogue with linear algebra, solutions to non-homogeneous non-linear systems are expressed through 

t:;j- ■ resultants by generalized Craemer rule [2j, also resultants control singularities of non-Gaussian integrals. 

For this reason, resultants have all chances to become the central special functions of emerging non-linear 
^^ ■ physics. They are already used in physical apphcations - see, for example, [3l|4l[5]. 

OO ! Practical application of non- linear algebra to scientific and engineering problems depends on development 

^D ■ of efficient methods to handle the resultants - as efficient as those for determinants in linear algebra. Like 

with determinants, not explicit expressions are needed (which take hundreds of pages already for low values 

of n and r's) but their properties and efficient computational algorithms, allowing to evaluate resultants in 

. , every particular situation. Several algorithms for calculation of resultants are already known [51 [71 [5] • In 

j^ ' this paper we describe yet another algorithm, which can be well computerized and used in practice. 

Our approach is based on the recursive relation 



logi?{/i, /2, ...,/„} = i . . . i log/i dlog/2 A ... A diogfn + ri logi?{/2, . . . , /„} 



xi=0 



(1) 



where the contour integral is (n — l)-fold, functions fi are taken with arguments /i(l, zi, . . . , z„-i) and the 
contour of integration encircles infinity. A sketchy proof of this relation is given in Appendix A. The relation 
of resultants to contour integrals of this kind is also studied in [7j, where promising connections with the 
theory of Riemann surfaces are pointed out. An explicit formula for \ogR follows after several iterations of 
([!]): for example, for n = 2 

l0gRr,rAflJ2}^ </log/i(l,z)dlog/2(l,z)-hrilog/2(0,l) 



for 



logRr^r2rAflJ2,f3} ^ f f log /l (1, ^1 , ^2)^ log /2 (1, ^1 , ^2) A d log /3(1, Zi , Z2) 

n /log/2(0,l,z)dlog/3(0,l,z) + rir2log/3(0,0,l) 



and so on. In this way, logarithm of the rcsuhant is expressed through the simple contour integrals. However, 
we would like to calculate resultants, not just their logarithms. To obtain R from log i?, we exponentiate 
the latter: for example, for n = 2 



Rr,r Ah J2}^expU log fl{l,z)d\ogf2{l,z)+n\ogf2{0,l) 



(2) 



Using this formula, it is straightforward to reproduce the well-known expression of the resultant of two 
polynomials through their roots: if 



/i(l,z) =aj|(z-aj) 
1=1 

h{l,z)^bl[{z-p,) 



then 

Rr 



(n 7-2 „ \ ri ?-2 

J2J2f log(^ - ^^)d^og{z - Pj) + r2 loga + n logfe = a'-^V' J] W^i " "*) 
i=li=l"' / i=lj = l 



Exponential formula ^ and its direct analogues for n > 2 have several advantages. The most important 
are simplicity and explicitness: such formulas have a transparent structure, which is easy to understand and 
memorize. Another one is universality: they can be written for any n, not just for low values. However, 
there is one serious drawback as well. Resultant is known to be a homogeneous polynomial of degree 



rir2 



di = 



in coefficients of the i-th equation /i, and of total degree 



rfi 



+ d„ 



rir2 . . . r„ 



1 1 

ri r2 



in coefficients of all equations, but in exponential formula the polynomial nature of the resultant is not ex- 
plicit. Therefore, it is not fully convenient for practical calculations of resultants (although it is conceptually 
adequate and can be used in theoretical considerations). 

We would like to convert the exponential formula into another form, which would be explicitly polynomial. 
Following [8], we suggest to apply H]) to a shifted system 

Mx) = ix2Y- - X2f2{x) 



. Jn\x) — \Xn) " — Xnjn\x) 

and expand the shifted logarithms in the integrands into Taylor series in powers of the "spectral parameters" 
Xi. We obtain the following series expansion 



logi?{/i, . . . ,/4 = - E • • • E ^fcifc2...fc„{/i, ■••,/«}• Af^A: 

fci=0 fc„=0 



.A^ 



(3) 



(a minus sign is just a convention) where particular Taylor components Tk-^k2...k„ are homogeneous polynomial 
expressions of degree ki in coefficients of ft. We call them traces 8] of a non-linear system /i(x), . . . , fn{x). 
Explicit formulas for traces are obtained in Appendix B, by making a series expansion in ([1]) and carrying 
out the remaining contour integrals. 

By exponentiating, we obtain the shifted resultant: 



i?{/i, ...,/„}= exp - ^ ... ^ Tfe,fe,...fc„ • A^ A^= . . . Af," 

\ fci=0 fc„=0 / 



The original resultant -R{/i, ...,/„} is equal to (—1)'' times the coefficient of Aj^^Aj^ 
We can extract it, expanding the right hand side in powers of A^: 



.Af,"ini?{/i,...,/„}. 



exP - E ■ • ■ E Tk,u,-^„ • At^ A^^ . . . A^ U ^ . . . ^ 7',,fc,...fc„ • A^A 

\ fci=0 fc„=0 / fci=0 fc„=o 



.A^ 



Such power series expansion of exp (5(x)), where S{x) is itself a power series (perhaps, of many variables), 
is often called a Schur expansion. It implies the following relation between Taylor components: 



fci+...+fe„ 



' k-ik^-.-kn 



E 



(-ir 



to! 



E 









^1 + ^2 + ■ ■ ■ + ^TTl — fc 



where the sum is taken over all ordered partitions of a a vector k — (fci, /c2, . . . , fc„) into m vectors, denoted 
as vi, . . . ^Vra- An ordered partition is a way of writing a vector with integer components as a sum of 
vectors with integer components, where the order of the items is significant. Polynomials V are often called 
multi-Schur polynomials. For example, 

Vifi = — Ti^o 

T'a.o = -T2,o + TlJ2 

^'2,1 = — T2,i + T2fiTo^i + Ti^oTis — 2^1,0^04/2 

^2,1,0 ~ ^T2s,o + T2fifiTos,o + 2^1,0,0^14,0 — 2^1, o,o2o, 1,0/2 

'^1,1,1 = ^2^1,1,1 + ri,o,o2o,i,i + 2o,i,ori,o,i + T'i,i,o2o,o,i ~ 2o,i,oTo,o,iTi,o,o 
Resultant is a certain multi-Schur polynomial of traces: 



R{fl,...Jn}^i-lfVd„ 



,d„ 



E 

m— 1 



(-1)^ 



to! 



E 






T,7„. 



Vi+V2 + ...+V^ 



(4) 



The explicit formula for traces, obtained in Appendix B, has a form 



2^ -Tk^k2...k„ 



fclfc2 






2<iJ<n 



(5) 



where 



(/)jij2,....j„ ~ coefficient of xj- 



Jl ^J2 



/(^ 



•2^n m /(,a:i,X2 



1 '^nj 



and the sum goes over all non-negative integer nx n matrices r^ such that the sum of entries in any i-th row 
or i-th coloumn is fixed and equals Viki. Formula ^ is valid for positive /ci, . . . , fc„. The cases when some 
ki — are equally simple, but they are somewhat degenerate and considered separately in the Appendix B. 



Together, formulas ^ and ([5]) constitute our final result, an explicit polynomial formula for the resultant 
of n homogeneous polynomials /i, ...,/„ in n variables. They provide a rather fast and effective computer 
algorithm for explicit evaluation of resultant for a finite set of polynomials fi . A MAPLE program, based 
on this algorithm, is presented in the appendix C. 



Appendix A: Proof of the relation ([T]) 

We begin from the well-known Poisson product formula [HIH]. It states, that the resultant of a system 

/i (a;i,a;2,. ..,x„) = 

f2{xi,X2, 



■ ; '^n } 

■ ; -^n ) 







(6) 



(7) 



fn (a;i,a;2,. ..,x„) = 
is equal to the product 

i?{/i,/2,...,/4 = C'.A(AW)A(A(2))...A(AW) 

where vectors A*-'-* are all the N ~ r2r^ ■ ■ -Tn common roots of polynomials /2, . . . , /«: 

' h (a(')) = 



^ In (a«) = 

This formula holds for the following reason. By definition of solvability, ^ is solvable iff /i vanishes on (at 
least) one vector A^'^ ^ 0. What is the same, resultant of ^ vanishes iff a product 

a(a(i))a(a(^))...a(aw) 

vanishes. Since resultant is a polynomial, it must be divisible by this product. The constant of proportionality 
C depends on the normalization of roots (but not on /i) and we assume that A*-*' are normalized in such a 
way, that C = 1. We proceed by rewriting this formula as 

R{h,h, ...,/„}= (A«Af ) . . . Al-')" . A (/i(i)) A (;i(^)) ... A (^-(-)) 

where /i^'^ — Iv^'^^ / h\' are the re-normalized roots, which have the first component jil' = 1. Since at the 
same time it follows from ([7]) that 



R{x^,h,...,fn}^h-rKY'...h-\ 



(1)a(2) 



AN) 



we have 



i?{/l, /2, . . . , /n} = R{xij2. . . . , InY' ' /l (/l^'^j /l (fi^^' 

Passing to logarithms, we obtain 

AT 



./i(m(^') 



logi?{/i,/2,...,/„} = rilogi?{xi,/2,...,/4 + El°g/i(/^^'0 



(8) 



In this form, Poisson product formula admits a contour integral representation. This is because the second 
term is a sum of values of one function at roots of the other (s). Using Cauchy residue theorem, such sums can 



be conveniently represented as contour integrals. Consider first an example, familiar from the elementary 
theory of functions of one complex variable. Let z be a complex variable, let f{z) and g{z) be polynomials 
of degrees |/| and \g\. Then, the sum of values of log/ at roots of g is given by 



^log/(piW) = j>logf{z)dlogg{z) 

i—l 



where the contour encircles all the \g\ roots of g, denoted by /i*^'-*, z = 1, 2, . . . , \g\. This is obvious, because 



\9\ 



ffW^ffoH^-^'^^ 



and 



i\ogf{z)d\ogg{z) ^Y^P^l^ dz = f]log/ (m(' 

i— 1 i—l 



This example is so simple, because any function of one complex variable is factorizable into linear factors. 
Such factorization has no analogue in the case of several complex variables, therefore, the proof is slightly 
more complicated in multidimensional situation. 

Proposition. The second term in the right hand side of ([8]) is given by a (rt — l)-fold contour integral 

^log/i (^^Wj == f ■■■ f log fi{l, Z2, ..., Zn)Yl dlog fkil, Z2, ■■■, Zn) 
1=1 •' ■' k=2 

where the contour encircles all the N = r2r^ ■ ■ -rn common roots of /2, . . . , /„, denoted by jl^'^'. 
Proof. Making a non-linear change of variables 

y2 = .f2(l,Z2,---,Zn) 
Vi = .f3{i,Z2,---,Zn) 



^ Vn = fnil, Z2, . . . , Zn) 



we obtain 



f ■■■ f logfi(l,Z2,...,Zn) Y[dlogfk{l,Z2,...,Zn) ^^ j> ■ ■ ■ ^ log /i \l,z^'^ [vYj d\ogy2 ■ ■ ■ dlogy^ 

k=2 i=l 

Since fk are polynomials of degrees r^, the inverse change of variables is multi- valued and has N = r2r3 . . . r„ 
branches, denoted as z^'^' (y) , i — 1,2, ... ,A^ in the above formula. The integral turns into a sum over 
branches. Integration over j/-variables is trivial 



log/i(l,z<^)(y)) 



dy2... dy„ 



log/i(l,z<^)(0)) 



y2---yn 

due to the Cauchy residue theorem. Obviously, the inverse change of variables maps (0,0,. . . ,0) into /i^^^: 

N N 



L.. /log/ill, Z2, . . . , z„) n d\ogh{l,Z2, . . . ,z„) = ^log/i (l,z<^) (o)) - ^log/i 

■^ '' k=2 2=1 i=l 



M 



(«) 



The proposition is proved. Therefore, we have 
logi?{/i,/2,...,/„} = rifogi?{a;i,/2,...,/„}+ f ■■■ f fog/i(l, ^2, • • • , z«) J]^ rflog/fc(l, 2:2, ■ • • , ^n) (9) 

•' •' fe=2 

Relation ^ follows from ^ and the simple identity 

R{xi,f2ixi,X2,. .. ,a;„),.. . , /„(a;i,X2, . .. ,a;„)} = i?{/2(0,2:2, . . . ,a;„), . . . , /„(0,a;2, . . .,a;„)} 
implied by the definition of resultant as a solvability condition. 



Appendix B: Explicit formula for the traces ([3] 

We are going to show, that for positive ki, the traces Tk-^k2...kn are given by 



Tk,k2...k„ == E , det (S.jrik, - r^) • J] ^^ 



where 



(/)iij2,...j„ = coefficient of xj'x^^ ■ • -a;^" in /(a;i,a::2, . . ■,Xn)'' 



for any polynomial /, and the sum goes over all non- negative integer n x n matrices rij such, that the sum 
of entries in any i-th row or i-th coloumn is fixed and equals riki'. 

n n 

J = l J = l 

Indeed, making in ([1]) a series expansion 

oo ^ 

log (zp - A,/,) = n logz. ^Y.T /^r''''A,* 
for positive ki we have 



k 

k=l 



Tk,M,...,kAf) = k.kZ'^l^ / ■ ■ • /^''"^ (/2^^2"'^^'^) A d (ftz^^'^) A ... A d {f!:-z-^-'-) 



Using 

2i + ... + iii— rifci 



/f*(l,Z2,...,Z„)= E (f^t^^ . Z^^.-.Z'^ 






we obtain after taking derivatives 

The Cauchy residue theorem implies a selection rule r^i + . . . + ri„ = r^fci. The transverse selection rule 



rii + . . .+r„j — Tiki is implied by homogenity of polynomials fi. Consider several examples in low dimensions: 



TkAD^E- 



fcl 



fcl 



Tki,k2,k3if) = J2 



Tk^MU) = E(^2fc2-r22) 

r2k2 - r22 -r23 

-r32 r^kz - r33 



fcl fc2 

(fl)''^ (f2)''" (fs)''" 

V^ -'-/riiri2r'i3 V-/ ^/r2ir22''23 v^'J/rait 

A;i fc2 A:3 



E 



^'2^2 - r22 -^23 -r24 

-'"32 ?'3^3 - ?'33 -7-34 

-'"42 -?'43 riki - r44 



J^ki,k2,k3,ki(j ) — 
V^i;riiri2ri3ri4 \j2jr21 



'•22l'23l'24 



l,/3Jr3ir32r33r34 W4Jr4i 



»"42l'43r44 



fcl 



fc3 



The cases, when some ki are equal to zero, also make no difficulty. If some ki = 0, then 

Tfci,fc2,...,fc„{/l,/2, ■ • ■ ,/n} = Ti Tk^,k2,...,k,-l,k, + l.---,k„{flif2, ■ ■ ■ ,/i-l,/i+l, ■ • ■ ,/n}|^.^o 

where the trace in the right hand side is taken in variables xi, . . . , Xi-i, Xi+i, . . . ,Xn- It is an easy exercise 
to prove this statement, making a series expansion in ([1]). 

Appendix C: MAPLE Program 



#################################### RESULTANT OF N HOMOGENEOUS POLYNOMIALS IN N VARIABLES 
interf ace(warnlevel = 0) : 
interf ace(displayprecision = 2) : 
with(combinat) : 
with(linalg) : 

##################### INITIAL DATA 

N := 3: 



f[l] 
f[2] 
f[3] 



= a * x[l]'2 + alpha * x[2]*x[3] 
= b * x[2]'2 + alpha * x[l]*x[3] 
= c * x[3]'2 + alpha * x[l]*x[2] 



################################### 

BeginTime := timeO: 

R := [seq(0, p = 1. .N)] : 

for i from 1 to N do 

if not ( selectCt -> op(0,t) = x, indets(f [i] ) ) subset {seq(x[p], p = 1..N)} ) then 
print( " Error: wrong variables. Must be ", seq(x[p], p = 1..N) ); 



return Error ; 
fi: 

R[i] := simplifyC ln( subs( {seq(x[p] = exp(l) * x [p] , p = 1..N)}, f [i] ) / f [i] ) ): 
if nops(indets(R[i] ) ) <> then 

print ( " Error: non-homogeneous polynomials " ); 

return Error ; 
fi: 
end do : 

deg := [seq( mul (R[pp] , pp = 1..N) / R[p] , p = 1..N)]: 
xDeg := add( deg[p], p = 1..N): 

Decompose := (u,v) -> {seq( [seq( composition(u+v,v) [p] [q] - 1, 

q = 1 . .nops(composition(u+v,v) [p] ) )] , p = 1 . .nops(composition(u+v,v) ))}: 

xDet := u -> if rowdim(u) = then return 1: else return det(u): fi: 

xTrace := proc(K): 

xl := select ( t -> K[t] = 0, [$1..N]): 
x2 := selectC t -> K[t] <> 0, [$1..N]): 
for p from 1 to nops(x2) do 

F[p] := - subs({seq(x[ss] =0, ss in xl) , x[x2[l]] = 1}, 
(f [x2 [p] ] ) -K [x2 [p] ] /x [x2 [p] ] - (R [x2 [p] ] *K [x2 [p] ] ) /K [x2 [p] ] ) : 
end do : 

dat := mul(R[ss] , ss in xl) * F[l] * xDet (matrix(nops(x2)-l ,nops(x2)-l , 
(i,j) -> diff( F[i+1], x[x2[j+l]] ) ) ); 
for p from 2 to nops(x2) do 

dat := eval(coeff( dat, x[x2[p]], -1)); 
end do : 

return simplif y(dat) : 
end proc: 

xxTrace := kk -> add( xTrace(vv), vv in Decompose(kk, N) ): 

xxSchur := kk -> add( mul( Tr[v[lll]], 111 = l..nops(v) )/nops(v)! * 

multinomial (nops (v) , seq( numboccur(v, convert (v, set) [i] ) , i = 1 . .nops (convert (v, set)))), 
V in partition(kk) ) : 

printC " Calculating resultant of type R[", op(R), "] " ); 
printC " Need to calculate ", xDeg, "traces " ); 
print ( " " ); 

for kkk from 1 to xDeg do 

Tr [kkk] := xxTrace(kkk) ; 

print ( kkk ) ; 
end do : 

printC " " ); 

printC " Resultant = ", factorC simplifyC C-l)~xDeg * xxSchur CxDeg) ))); 
printC " Done in ", timeC) - BeginTime, " sec. "); 
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